Density gradients for the exchange energy of electrons in two dimensions 
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We derive a generalized gradient approximation to the exchange energy to be used in density 
functional theory calculations of two-dimensional systems. This class of approximations has a long 
and successful history, but it has not yet been fully investigated for electrons in two dimensions. We 
follow the approach originally proposed by Becke for three-dimensional systems [Int. J. Quantum 
Chem. 23, 1915 (1983), J. Chem. Phys. 85, 7184 (1986)]. The resulting functional depends on two 
parameters that are adjusted to a test set of parabolically confined quantum dots. Our exchange 
functional is then tested on a variety of systems with promising results, reducing the error in the 
exchange energy by a factor of 4 with respect to the simple local density approximation. 
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I. INTRODUCTION 

Present nanoscale electronic devices contain a large va- 
riety of low-dimensional systems in which the many-body 
problems of interacting electrons need to be addressed. 
These systems include, e.g., modulated semiconductor 
layers and surfaces, quantum Hall systems, spintronic de- 
vices, and quantum dots.^ In order to describe the elec- 
tronic properties of these systems, a practical and accu- 
rate way of computing the energy components is required. 
Since the advent of density-functional theory^'^''* (DFT) 
much effort went into the development of approximate 
functionals for the exchange and correlation energies. 
Most of this work focused on three-dimensional (3D) sys- 
tems, where considerable advances beyond the commonly 
used local density approximation (LDA) were achieved 
by generalized gradient approximations (GGAs), orbital 
functionals, and hybrid functionals.^ However, previous 
studies have shown that most functionals developed for 
3D systems break down when applied to two-dimensional 
(2D) systems 

Within the DFT approach, 2D systems are most com- 
monly treated using the 2D-LDA exchange [see Eq. (13) 
m Ref^, which is then combined with the 2D-LDA corre- 
lation parametrized first by Tanatar and Ceperley^ and 
later, for the complete range of coUinear spin polariza- 
tion, by Attaccalite and co-workers>i^ Despite the rel- 
atively good performance of LDA with respect to, e.g., 
quantum Monte Carlo calculations,^^ there is a clear lack 
of accurate 2D density functionals. 

The exact-exchange functional employed within the 
optimized effective potential method, which automati- 
cally conforms to various dimensionalities, seems an ap- 
pealing alternative to the LDA, and it has recently been 
applied to quantum dotSjiS. In that method, however, the 
development of approximations for the correlation ener- 
gies compatible with exact-exchange energies remains a 
complicated problem. 



The derivation of local and semi-local approximations 
for the exchange-correlation energy in 2D can be car- 
ried out following the general lines already employed for 
the 3D case. In this way one can take advantage of the 
almost 40 years of experience in that field. Only very 
recently, such efforts have been done in 2D by develop- 
ing exchange functionals specially tailored for finite 2D 
systemsi ^3ii4, 15^16 ^^lis work we take the most natu- 
ral step beyond the LDA by including the density gra- 
dients in the functional. To this end, there are several 
possible approaches, e.g., the gradient expansion of the 
exchange holer^ semiclassical expansions from the Dirac 
or Bloch density matrix;^ and the GGAs of Perdew^^ and 
Beckeji2i2£iiSi Here we follow the approach introduced by 
Becke for 3D systems in Refs.[il| and [2^ and derive and 
apply a GGA for the exchange energy of 2D electronic 
systems. Tests for a diverse set of 2D quantum dots show 
excellent performance of the derived approximation when 
compared with exact-exchange results. 



II. DERIVATION OF THE APPROXIMATIONS 

Within the Kohn-Sham approach to spin-DFT,^- the 
ground state energies and spin densities ^^(r) of a sys- 
tem of iV — N-^ + Ni interacting electrons are deter- 
mined. The total energy, which is minimized to obtain 
the ground-state energy, is written as a functional of the 
spin densities (in Hartree atomic units) 



dVw^(r)n^(r), (1) 



where T's['T'cr] is the Kohn-Sham kinetic energy functional, 
Va- (r) is the external (local) spin-dependent scalar poten- 
tial acting upon the interacting system. En [n] is the clas- 
sical electrostatic or Hartree energy of the total charge 
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density n(r) = n|(r)+n|(r), and i?xc[«cr] is the exchange- 
correlation energy functional. The latter can be further 
decomposed into the exchange and correlation parts as 

In this work we focus in the exchange-energy func- 
tional, that can be expressed as 



. (2) 

where, within the restriction that the nonintcracting 
ground state is nondegenerate and hence takes the form 
of a single Slater determinant, the exchange- hole (or 
Fermi-hole) function is given by 

= ^) ■ 

The sum in the numerator is the one-body spin-density 
matrix of the Slater determinant constructed from the 
Kohn-Sham orbitals, fk.a- Moreover, integrating this 
function over r2 yields 



Jd^r2 /ix,^(ri,r2) = 1. 



(4) 



This exact property reflects the fact that around an elec- 
tron with spin a at ri, other electrons of the same spin 
are less likely to be found. This a consequence of the 
Pauli principle. From Eq. ^ it is clear that to evalu- 
ate the exchange energy in 2D, we just need to know the 
cylindrical average w.r.t. s = r2 — ri of the exchange-hole 
around ri, 
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/ix,CT(ri,ri +s). 



from which 
E^[na] - 



-TT E ld^rna{r) /ds /ix,<T(r, s). 



(5) 



(6) 



where we have renamed ri as r. Expressing the exchange 
hole by its Taylor expansion, and considering its cylin- 
drical average, one arrives at the following expression. 



(7) 



where is the so-called local curvature of the exchange 
hole around the given reference point r. This function 
can be expressed a o^'^'^^i^^'^'^ 



where 



V^n^ir) - 2T,(r) 



1 |V7i,(r)r 



(8) 



(9) 



is (twice) the spin-dependent kinetic-energy density, and 



ip,4r) = ^ £ {^fe,.(r)[V(^fe,.(r)]- 

k=l 

[VvlAr)]^kA^)} (10) 

is the spin-dependent paramagnetic current density. 
Both To- and jp.o- depend explicitly on the Kohn-Sham 
orbitals. Thus the expression in Eq. ([5]) has an implicit 
dependence on the spin densities n^. 

A. Small density-gradient limit 

When the inhomogeneity of the electron system is 
small, we may regard the homogeneous 2D electron gas 
(2DEG) as a good reference system. In this case, we have 
an exact expressio: 



,25 



■;2DEG 



(11) 



where kp^^r = V??™^ is the Fermi momentum (for spin 
a) in 2D, and Ji is the ordinary Bessel hmction of the 
first kind in the first order. Notice that the principal 
maximum of Eq. pT|) accounts for 95% of the exchange 
energy. We thus follow the idea introduced by Becke for 
the 3D case^ to modify the principal maximum of h^J^'^ 
by a polynomial factor. This improves the short-range 
behavior but leaves the secondary maximum unchanged. 
We may write 

/ix,.(r, s)=[l + a,{v)s^ + 6,(r)/ + . . .] hf^'^is), 

(12) 

for kp^^s < z, and 



(13) 



for kp^cS > z, where z is the first zero of Ji. 

Now we have to find an expression for both and 
ba- Comparing Eq. with Eq. ([7]), making use of the 
series representation of Ji 



^i(y)-E 



k=0 



(-l)'(f)' 

A:!r(fc + 2) 



(14) 



and replacing the expression of t„ in Eq. Q by the 2D 
Thomas-Fermi expressio n^'^i^^i^^ 



r.(r)=27rn^(r) + iW(r) + %^, 
3 ncr(r) 



we arrive at 



k=l 



4n^(r) 



(15) 



(16) 
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Next, we can determine ba by applying the normalization 
constraint of Eq. ^ , leading to 



/(3) 

Here the symbol I{m) denotes 



I(m) 



dyy"Jf(y), 



(17) 



(18) 



where z is again the first zero in Ji. The values of the 
integrals I{n) can be determined numerically. 
Making use of Eq. ^ , we arrive at 

[n,] = £;2DEG ^ ^SGL (19) 

where 



with 



1 



43/2, 



/(0)/(3)-J(l)/(2) 
/(3) 



(20) 



(21) 



where SGL refers to the small-gradient limit. Using 
Green's first identity when integrating the first term, we 
find 



E: 



SGL 



«SGL h?rnl'\v)xl{r). 



with KsGL = 5k/6 and Xc^ir) = |VnCT(r)| /'n^J'^{v) is the 
usual dimensionless parameter for exchange. 



B. Large density-gradient limit 

Using Eqs. ^ and (fT5|) . we can rewrite Eq. ([7]) as 



1 |Vrv(r) 

2 n^(r) 



(23) 



Here we follow the reasoning of Becke^S (applied in 3D), 
and assume that the term in the density gradient domi- 
nates over the other terms. Thus, we obtain 



h (r.)«i^^.^ 



(24) 



We note that expression ([M]) is valid only for small s. 
Otherwise, following again the argument of Becke, we 
propose 



irvrv(rr^2 

8 n,{v) 



F{a,{r)s), (25) 



where 



Fiy) 



(26) 



The form of F in Eq. (pS)) corresponds to a Gaussian ap- 
proximation for the exchange-hole. This would be exact 
in the case of a single electron in a harmonic confinement 
potential. However, another choice for F may be consid- 
ered; so far it allows to reproduce the correct short-range 
behavior given by Eq. (j24p . and decay in a way leading 
to finite exchange energies. 

In Eq. (j25p . the parameter can be determined by 
enforcing the normalization condition of Eq. (jj]). This 
leads to 



TT G(3) |Vn^(r)r 
4 n^(r) 



where 



G(m) ^ / dyy'^e-y . 
Jo 

Finally, making use of Eq. we arrive at 



(27) 



(28) 



where 



23/2 



d2rn3/2(r)4/2(r), (29) 



(30) 



G(2)G-3/4(3) 



(22) Here LGL refers to the large-gradient limit. 



III. EXCHANGE-ENERGY FUNCTIONAL 

We can combine Eqs. (|22|) and ([29]) in an expression 
that interpolates the exchange energy of an inhomoge- 
neous electron gas from the small density-gradient limit 
to the large one. A possible expression is of the form 



,3/2 



(r) 



[l + jxl{r)f^' 



(31) 



or, writing this expression explicitly in terms of the den- 
sity and its gradients. 



E^'^^[n,,Vn,] = El^^'^in,] 



|Vn„(r) 



nj (r) 



1+7^ 



(r) 



3/4 



(32) 



where (5 and 7 can be trivially written in terms of klgl 
and KsGL- At last, to improve the flexibility of the above 
approximate functional, we replace both (3 and 7 with 
two parameters to be determined by fitting the exchange 
energies of an ensemble of physically relevant 2D systems. 
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FIG. 1: (color online) Exchange energy per electron for a 
series of parabolically confined quantum dots with A'' electrons 
and confinement strength u. For clarity, the line for N — 6 
was shifted down by 0.5 a.u., the line for = 12 by 1 a.u., and 
the line for A'^ = 20 by 1.5 a.u. The lines are: exact exchange 
(solid gray), exchange-only LDA (dotted red), and this work 
(dashed green). 



TABLE I: Exchange energies (in atomic units) for parabolic 
quantum dots. The last row contains the mean percentage 
error, A. The first four lines represent the systems that were 
used in the fitting of l3 and 7. The columns are: exact ex- 
change (EXX), local density approximation (LDA), and the 
generalized gradient approximation presented in this work 
(GGA). 



N 






TTiLDA 
--fix 




2 


1 


1.083 


0.9672 


1.051 


2 


1/4 


0.4850 


0.4312 


0.4704 


2 


1/16 


0.2073 


0.1843 


0.2023 


2 


1/36 


0.1239 


0.1108 


0.1276 


6 


0.42168 


2.229 


2.110 


2.206 


6 


1/1.89^ 


1.735 


1.642 


1.719 


6 


1/4 


1.618 


1.531 


1.603 


12 1/1.89^ 


3.791 


3.668 


3.777 


A 






7.9% 


1.8% 



TABLE IL Exchange energies (in atomic units) for square 
(area = tt x tt) quantum dots. The meaning of the columns 
is identical to Table |T] 



N 


^EXX 


771LDA 


_^GGA 


2 


1.417 


1.288 


1.383 


6 


6.147 


5.902 


6.180 


8 


9.509 


9.017 


9.434 


12 


16.24 


15.91 


16.46 


16 


25.23 


24.35 


25.15 


A 




4.8% 


1.1% 



N in the dot and the confinement strength uj. From the 
fully self-consistent DFT results summarized in Fig. [1] it 
is clear that our functional outperforms the LDA for the 
whole range of N and w. A more quantitative picture 
can be obtained from Tables [J and [IT] where we present 
numerical values for the exchange energy obtained with 
different approximations, for both parabolically confined 
and hard- wall square^ quantum dots, respectively. It is 
clear that our functional yields errors that are smaller by 
at least a factor of 4 than the errors of the simple LDA. 

We tested the performance of our GGA functional also 
in the large- iV limit. Here we used the exact exchange 
energies known for parabolic closed-shell quantum dots 
of noninteracting electronsi^ Using a fixed, analytic elec- 
tron density as an input in Eq. (j32p . we found relative 
percentage errors of 0.7% and 0.5% for N = 110 and 
420, respectively, whereas the corresponding errors of 
the LDA are 0.5 % and 0.2 %. Even if the plain LDA 
is slightly more accurate in this case, we find the result 
very satisfactory in view of the fitting of (3 and 7 accord- 
ing to two-electron data as explained above. We point 
out that in terms of potential applications for our GGA 
functional, the main interest is in relatively small sys- 
tems (up to a few dozen of electrons). However, the 
GGA described here - possibly with simple modifica- 
tions - enables precise many-electron calculations also 
in large-scale structures such as quantum-dot arrays^! or 
quantum-Hall devicesi^ 



We have chosen to fit our parameters to a set of four 
parabolic two-electron quantum dotai with confinement 
stengths w = 1, 1/4, 1/16, and l/36a.u., respectively. 
These are very well studied systems^^ that span a wide 
range of the density parameter, 0.9 < < 5.7, where 
Ts = N~-^^^ uj~^^^ . For the reference exchange energies 
we performed self-consistent exact-exchange calculations 
with the code octopus^ using the Krieger-Li-Iafrate^ 
(KLI) approximation, which is a very accurate approxi- 
mation in the static case.'^'- We obtained /? = 0.003317 
and 7 — 0.008323, which are close to the parameters 
found by Becke for the 3D case by fitting a series of rare- 
gas atoms, = 0.00375 and 7^0 = 0.007 (Ref. [13). 

Our functional was then applied to parabolically con- 
fined quantum dots by varying the number of electrons 



IV. CONCLUSIONS 

We proposed a generalized gradient approximation for 
the exchange energy of two-dimensional systems, follow- 
ing the same lines developed for three-dimensional sys- 
tems proposed by Becke . ^^1^° Analyzing the small- and 
large-density gradient limits, we arrived at the expression 
of a functional which depends on two parameters, that 
are fitted to a test set composed of four parabolically 
confined quantum dots containing two electrons. Further 
calculations, both for parabolically confined and square 
quantum dots, show that our approximations yields er- 
rors that are at least a factor of 4 better than the local 
density approximation. 
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We believe that this is the necessary step in the con- 
struction of rehable generahzed gradient approximations 
for two-dimensional systems. The next steps will neces- 
sarily involve the construction of correlation functionals 
(beyond the local density app roximation^" or the forms 
proposed in Refs. [l^ and Il6l). improvements to the ex- 
change functional presented in this work, and further ex- 
tensive tests, and analysis, to assess the quality of the 
same functionals. In this path we expect that the expe- 
rience in the development of functionals in three dimen- 
sions will be a very useful guide, but we can also expect 
that experience gained in describing exchange and cor- 
relation in these low-dimensional systems can be again 
transposed to bring new insights and ideas into the three- 



dimensional world. 
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